Coexistence of Anomalous and Normal Diffusion in Integrable Mott Insulators 
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We study the finite-momentum spin dynamics in the one-dimensional XXZ spin chain within the 
Ising-type regime at high temperatures using density autocorrelations within hnear response theory 
and real-time propagation of nonequilibrium densities. While for the nonintegrable model results 
are well consistent with normal diffusion, the finite-size integrable model unveils the coexistence of 
anomalous and normal diffusion in different regimes of time. In particular, numerical results show 
a Gaussian relaxation at smallest nonzero momenta which we relate to nonzero stiffness in a grand 
canonical ensemble. For larger but still small momenta normal-like diffusion is recovered. Similar 
results for the model of impenetrable particles also help to resolve rather conflicting conclusions on 
transport in integrable Mott insulators. 
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Introduction. — Theoretical investigations of transport 
in many-body systems of interacting fermions estab- 
lished several novel, entirely quantum aspects which go 
well beyond usual weak-scattering or Boltzmann-type ap- 
proaches to transport. It has been shown that integra- 
bility of a model system can change qualitatively the re- 
sponse to external driving. A prominent experimentally 
relevant example of such a system is the one-dimensional 
(ID) Heisenberg (XXZ) model Q, where a long-standing 
question is the existence of diffusion in the isotropic case 
[2-|Jl , and recently also spin systems mapping on the 
anisotropic (Ising-like) case became of interest [5[. 

While it is by now quite well established that inte- 
grable conductors, in the easy plane regime, exhibit at 
any temperature T > ballistic (dissipationless) trans- 
port within linear response (LR) theory 0] charac- 
terized by a finite spin stiffness Ds{T > 0) > 
which was recently confirmed by exact lower bounds 
T > transport in the Ising-type (Mott insulating in the 
fermionic representation) regime of the same model still 
represents a challenge with some apparently confiicting 
conclusions. While Ds{T > 0) ^ ioT an infinite system 
in this regime, numerical studies for the LR dynamical 
conductivity cr^uj) at high T ^ reveal in finite sys- 
tems of length L a broad incoherent response with on 
the one hand quite featureless spectra in the thermody- 
namic limit L — ?> oo and consequently a finite d.c. value 
CTdc = f^i^ 0) [ll]. On the other hand, the low- 
Lu behavior is dominated by a huge finite-size anomaly 
with vanishing response within a window cj < w* oc l/L. 
Anomalous behavior reappears also at finite-field driv- 
ing or weak perturbation (of the model), which both 
break the integrability and indicate on the existence of an 
"ideal insulator" with adc — >■ [14| as the proper limit. 
Also the exact nonequilibrium steady state of a strongly 



driven open XXZ chain [12| reveals a similar anomaly 
with the current decaying exponentially with the length 
of the chain. This is in contrast with steady transport 
under near- equilibrium conditions suggesting again a fi- 
nite diffusion constant V which is also consistent 
with previous studies performed at finite time [l6| and 
momentum [17j |. 

The aim of this Letter is to reconcile apparently in- 
consistent manifestations of diffusion in the (Mott) in- 
sulating regime at T > 0, in particular at high T — oo, 
whereby we concentrate our analysis on the XXZ Heisen- 
berg model with the anisotropy A > 1. We first note 
that anomalous transport in the integrable model can 
be related to the absence of a characteristic scale repre- 
senting "the mean free path" I*, which is in small sys- 
tems substituted effectively by the actual size L. On the 
other hand the diffusion constant V in the thermody- 
namic limit L — )■ CX3 would indicate a very short I* ^ 1. 
This dichotomy shows up for finite systems in spin corre- 
lations Sq{t) at momentum g > as diffusion- type decay 
at t < t* while at t > t* the decay becomes Gaussian, 
whereby t* oc L. The latter behavior is shown to be dom- 
inant at smallest non-zero q = qi = 2tt/L and its origin 
can be traced back to the existence of a finite stiffness 
Ds > in a grand canonical ensemble. With increas- 
ing q > qi normal diffusion prevails. To strengthen our 
arguments and results we confirm the same phenomena 
within the ID model of impenetrable particles (U ^ oo 
Hubbard model), where clearly no steady spin current 
is possible at zero magnetization, nevertheless S'g>o(t) 
again reveals a coexistence of normal diffusion and Gaus- 
sian decay. 

ID Heisenberg model. — First, we address the question 
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of spin transport in the ID anisotropic Heiscnberg model, 

L 
r=l 

(1) 

where 5* {i = x,y,z) are spin s = 1/2 operators at 
site r, L the length of the chain with periodic bound- 
ary conditions (p.b.c), and A represents the anisotropy. 
We allow also for a next-nearest neighbor z2-interaction 
with A2 7^ breaking the integrability of the model. 
It should be reminded that the Hamiltonian ([T|) can 
be mapped on a t-V-W model of interacting spinless 
fermions with hopping t ~ J/2 and inter-site interac- 
tions V = J A, W = JA2. In this fermionic picture, spin 
transport corresponds to charge transport and the here 
interesting regime of A > 1, S'^ot = (at A2 ^ 0) to the 
Mott insulator. 

Our aim is to analyze spin transport in the insulating 
regime A > 1 at T 0. Two complementary numerical 
approaches are used: a) LR theory calculating relevant 
dynamical spin correlation functions and b) real-time 
propagation (TP) of spin after switching off a perturbing 
magnetic field. The latter can be regarded as the con- 
trol of the validity of LR theory at finite perturbations, a 
question being nontrivial in particular for integrable (and 
nonergodic) systems. In both approaches, however, we 
apply essentially the same numerical approach by means 
of Lanczos diagonalization of finite systems [Til l , applica- 
ble up to L ~ 30, beyond the range of full diagonalization 
(FD), where L<20^^. 

Within the framework of LR theory we consider 
the time-dependent spin correlation function Sq{t) = 
Re(S'|(t)5iq)/L where = E^e"''''?^ and, due to 
p.b.c, q = 2'Kk/L. Here, (...) denotes the thermody- 
namic average at temperature T. We primarily focus 
on the high-temperature limit /? = 1/T — > at zero 
magnetization S^^^ — where we use the microcanoni- 
cal Lanczos method (MCLM) [l^l to evaluate ^^(a;) in 
finite systems and then perform the Fourier transform 
into t-dependent Sq{t). In the case of perfectly diffu- 
sive dynamics we would expect Sq(t) oc exp(— g^I^i). 
In general, however, the instantaneous rate 'Dq{t) = 
~ S q(t) / [q^ S q{t)] can become constant only in a hy- 
drodynamic regime at small enough q and long t. For 
such q, 2?<j(i) is related to the autocorrelation func- 
tion Jo(t) = (J§{t)J§}/L of the q = current Jn = 
JT,riSrSr+i - 3^+^) by the Einstein relation (at 
/3^0) 

MmVqit) = ^ = ^ fdt' Mt') , (2) 
9^" X So(t = 0) Jo 

assuming non-singular behavior at g — > 0. This assump- 
tion becomes relevant since a finite system features a non- 
vanishing stiffness Ds, i.e., Jo{t > t*) ~ 2Ds, in partic- 
ular if one considers the grand canonical averaging (over 



all in the XXZ model) where D,, oc l/L as dis- 
cussed furt heron. Hence, due to 

lim P,(t>r)= const. + ^1^, (3) 

a finite stiffness Dg restricts normal diffusion to t < t* 
and implies anomalous dynamics for t > t* [lij . 

To study spin transport using the real-time dynam- 
ics, we extend the Hamiltonian ([l} by introducing a 
position- and time-dependent magnetic field, H ^ H — 
X^r hr{t)S^. The initial equilibrium state corresponding 
to small but finite inverse temperature /3 -C 1/ J is then 
obtained by means of the MCLM [l^ with the initial per- 
turbation hr{t < 0) = ft.ocos(gr). Finite magnetic field 
induces site-dependent magnetization which for /3 <C l/J 
is on average {S^){0) « /3/ir(0)/4. At t = the field 
is switched off and the system evolution is obtained 
from the Lanczos time propagation method [13, HM . 
Following the relaxation of local magnetization {S^){t) 
we first check that the magnetization preserves its ini- 
tial spatial profile for the assumed hr{t), i.e., {S^){t) = 
(5^)(0)/(i) even when the dynamics are anomalous and 
strongly deviate from a simple exponential diffusion de- 
pendence {S^){t) = {S^){0)c^p{~VqH). Therefore, f{t) 
can be used for distinguishing between normal diffu- 
sion and anomalous dynamics. For /lo — >■ this quan- 
tity should be compared with the ratio Sq{t)/ Sq{0) from 
the LR approach. Such a comparison for smallest finite 
q = qi is shown in Fig. [T] A quantitative agreement 
between both methods is clearly visible, confirming that 
the LR approach (for q > 0) remains for finite valid 
even when the system is integrable. On the other hand, 
the key result in Fig. [1] concerns clear presence of normal 
diffusion in a generic nonintegrable system (A2 ^ 0) in 
sharp contrast with fast and anomalous relaxation visible 
in the integrable system. 

Next we focus on the behavior of the integrable case 
with A2 = 0. To give insight into the origin of the fast 
decay of Sq(t) we consider the instantaneous rate 'Dq{t). 
In Fig. [2] we thus summarize our numerical LR results 
on 'Dq{t) at q = qi and A = 2.0 for chains ranging from 
L = 14 to 30, where FD is used for i < 20 and the 
MCLM for L = 24-30. Apparently, Vq{t) first increases 
at short times t ^ 1.5/ J and then develops a rather con- 
stant plateau r'g(i) ~ 0.4 J at intermediate times t < t*, 
which is consistent with previous studies at finite momen- 
tum yielding Vq{t) « 0.88J/A Clearly, the plateau 
marks diffusive dynamics at intermediate time scales and 
we further observe this time scale to increase with sys- 
tem size approximately as t* w L/{3J), see Fig. [5] (inset). 
While this scaling with L is a pointer to purely diffusive 
dynamics in the thermodynamic limit L — > 00, the dy- 
namics at long times t > t* turns out to be different 
for finite L, in particular for q = qi. As clearly visible 
in Fig. [5J 'Dq{t) increases linearly with time, which in- 
dicates anomalous dynamics. One might be tempted to 
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FIG. 1. (color online) Decay of spin densities Sq{t) in the 
anisotropic Heisenberg model with A = 1.5 shown for smallest 
nonzero q — qi on a. chain of length L = 26 and A2 = 0, 
0.5, corresponding to integrable and nonintegrable models, 
respectively. LR results are obtained via the MCLM method 
at /3 = while TP results are obtained at /3 = 0.4/ J. 



FIG. 3. (color online) Decay of densities Sq vs. scaled time 
tq^ at anisotropy A = 2.0, calculated within LR via MCLM 
on chains with L = 24, 28 sites for different wave vectors 
q — qi, 2qi , respectively. Dashed curve indicates exponential 
decay with D = 0.4J. The position of the marginal time t* 
(= 2) is also marked for L = 28. 



relate the linear increase directly to a finite stiffness Ds, 
cf. Eq. (jSj. But in a canonical ensemble at zero mag- 
netization the stiffness decreases exponentially fast with 
system size [l3| while the slope in Fig. [2] scales rather 
as 1/L. Hence, we compare the slope with the stiffness 
resulting for a grand canonical ensemble (over all S^^^ 
sectors) with zero average magnetization. Since lower 
bounds for Dg are given by the Mazur inequality, and 
a projection to the conserved energy current has been 
shown to represent well the actual stiffness we ar- 
rive at 



A2 J2 



4(1 + 2A2)L' 



(4) 



and, noting the sum rule Salt = 0) = 1/4, we obtain 
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FIG. 2. (color online) Spin decay rate T>q{t) at anisotropy 
A — 2.0, A2 — calculated for the smallest finite gi in chains 
of different L — 14, 20, 24, 26, 30. Dashed curves indicate 
straight lines oc l.OtJ^/L while t* = L/{3J) is marked for 
L = 30. The inset emphasizes the L-dependence at short 
times. 



2Df/Sa{t = 0) > 0.89 JVL for A = 2.0 while FD results 
yield « I.IOJ^/L for small system sizes 0, [H]. The 
convincing agreement with the slope of 'Dq{t) in Fig. [21 
K, 1.0J2/L, is a hint at an effective finite stiffness at 
nonzero q — qi, not being reported yet. In any case, 
the linear increase of the rate Vq (t) at longer times t > 
t* identifies the fast Gaussian decay of Sq{t) in Fig. [U 
similarly as found for transport in complex one-particle 
models of finite size 
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Even though not shown here 
explicitly, this type of relaxation also manifests in the 
spectrum Sqicj) as an anomaly of Gaussian shape at low 
frequencies oj < cj* ~ 27r/t* oc 1/L, well pronounced at 
q = qi since in this case the main part of the sum rule is 
located at ut < ui* . 

So far, we have concentrated on Sq{t) at the small- 
est finite q = qi. Here, the effective relaxation time 
Tq ^ t* and the exponential decay only appears as a 
minor fraction of the total relaxation while Gaussian re- 
laxation dominates, see Fig. [3] On the other hand, for 
large enough g > gi, we realize that t ^ t* and the ex- 
ponential relaxation starts to dominate the decay curve. 
Thus, a pertinent criterion for "normal" diffusion relax- 
ation is given by r = l/{q^T>) < t* . The latter criterion 
already turns out to be rather well satisfied for q = 2qi 
for the considered chain lengths, e.g., for L = 28, as il- 
lustrated in Fig. [3] 

ID model of impenetrable particles. — We are finally 
going to address the question to which extent the ob- 
served dynamical behavior are generic for integrable 
quantum systems in the insulating (Mott-Hubbard-type) 
regime. To this end we investigate the ID model of im- 
penetrable particles, which has also been shown to be- 
have anomalously with respect to transport [l3| . The 
model (being the U ^ 00 limit of the Hubbard model or 
J ^ limit of the t-J model) is given by the Hamilto- 
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nian 

L 

^ = -^~EE'^'+i,A.+H.c., (5) 

r— 1 s 

where projected fermion operators Cr,s = Cr,s(l — ?^r,-s) 
take into account that double occupancy of sites is for- 
bidden. The two different species of particles arc given 
by up (t) and down (X) spin fermions. Wc should note 
that there is a close analogy of the XXZ model ([1]) in the 
large anisotropy (Ising) limit A ;2> 1 with the i model 
([S]). Namely, within the Ising limit we are dealing with 
the Neel ordered ground state (at S^^^ = 0) and the 
excited states composed of split subspaces of oppositely 
charged "soliton-antisoliton" (ss) pairs. In such a limit, 
the solitons/antisolitons behave effectively as impenetra- 
ble quantum particles since their crossing would require 
virtual processes with energy 5E ~ JA within the XXZ 
model (or 5E ~ U within the Hubbard model). 

Within the t model, the charge and spin currents can 
be written as 

L 

jlc.s] ^ s] 4+i,.5... + H.c. (6) 

r— 1 s 

with = [2,l/2]n(l - where n = 

(N^+Ni) / L is the filling. Further we notice that J§ com- 
mutes with H (from the perspective of charge the model 
is equivalent to ID noninteracting spinlcss fermions) 
while Jq docs not. Their overlap {JqJq)/L — 2to(1 — n)P 
vanishes when magnetization m = (N^ — N\^)/{2L) is 
zero. Moreover, it is quite evident that for to = there 
could be no d.c. spin transport since particles can- 
not cross with particles which implies Dg = but 
as well CTdc = cr(w — > 0) = [l^. On the other hand, 
using a grand-canonical ensemble with average to = 0, 
the Mazur inequality leads to the lower bound 



Noting the sum rule 5*0 (t = 0) = n/4, we obtain 
2Df/So{t = 0) > 2.0P/L for n ^ 1/2 (corresponding to 
quarter filling for the Hubbard model), as considered in 
the following. This lower bound we again compare with 
the instantaneous rate 'Dq{t) at smallest q ~ qi in a chain 
of length L = 20, maximally treatable with MCLM. As 
shown in Fig. 2] (inset), we indeed find 'Dq{t) to increase 
remarkably well linearly at long times t > t* , as before, 
in obvious agreement with the lower bound. This linear 
increase leads to a Gaussian decay of Sq {t) at long times 
t > t* . On the other hand, frequency moments and the 
limit i — >■ oo of the dynamical spin conductivity (t(w) are 
again consistent with a finite diffusion constant being in 
this (n = 1/2) case V = 0.76i [l^. Also we observe the 
dynamics at t < t* to be consistent with an exponential 




FIG. 4. (color online) Decay of spin density autocorrelations 
Sq{t) in the t model at n = 1/2 (quarter filling), zero magne- 
tization, and high temperatures, depicted for the two smallest 
momenta g > in a chain of length L = 20. One solid curve 
indicates an exponential decay with a constant rate U = 0.76t. 
Inset; Rate T>q(t) at the smaller of both q for the same pa- 
rameter set. The dashed curve indicates a line oc Q.ltP . 



decay involving this value for the diffusion coefficient, as 
shown in Fig. 2] As well we confirm in Fig. 2] that for 
larger q = 2qi the decay approaches the "normal" diffu- 
sion behavior. 

Conclusion. — In summary we studied the finite-g spin 
dynamics in the ID Hcisenberg chain with anisotropy 
A > 1 in the high-temperature limit /? — 5- 0. As one of 
the main results, we first showed the validity of linear re- 
sponse theory at finite perturbations using the real-time 
propagation of nonequilibrium densities. While wc found 
exponential relaxation (normal diffusion) in the noninte- 
grable model, we observed in the integrable model the co- 
existence of a Gaussian relaxation (anomalous diffusion) 
at long times t > t* oc L, being dominant at smallest 
q ~ qi where the effective relaxation time of spin modu- 
lations is Tq ^ t* . On the other hand, when increasing 
q > qi, normal diffusion prevailed and also the respec- 
tive diffusion constant is in quantitative agreement with 
transport coefficients from steady state scenarios [l3| ■ To 
be in full agreement with the latter (open system) sce- 
narios it is therefore irnportant to perform the limits in 
the appropriate order [22l |. i.e., first L oo {t* ^ oo) 
and then q ^ (r^ — oo), although the opposite limits 
can be as well relevant and realized [31 ■ Finally, we ob- 
tained similar results on the ID model of impenetrable 
particles, suggesting that the observed dynamics is quite 
generic for integrable Mott insulators. 
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